Ground State Fidelity in One Dimensional Gapless Models 
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A general relation between quantum phase transitions and the second derivative of the fidelity (or 
the "fidelity susceptibility") is proposed. The validity and the limitation of the fidelity susceptibility 
in characterizing quantum phase transitions is thus established. Moreover, based on the bosonization 
method, general formulas of the fidelity and the fidelity susceptibility are obtained for a class of 
one-dimensional gapless systems known as the Tomonaga-Luttinger liquid. Applying these formulas 
to the one-dimensional spin-1/2 XX Z model, we find that quantum phase transitions, even of the 
Beresinskii-Kosterlitz-Thouless type, can be signaled by the fidelity susceptibility. 
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In order to obtain fresh insight into the quantum 
many-body problem, a great deal of effort has been de- 
voted to the application of the concepts emerged from 
quantum information science^ to the analysis of quan- 
tum phase transitions (QPTs).— QPTs are characterized 
by the drastic change in the ground states of quantum 
many-body systems, driven solely by quantum fluctua- 
tions. Since entanglement measures the strength of quan- 
tum correlations between subsystems of a compound sys- 
tem, it is natural to expect that entanglement will be 
a reliable indicator of QPTs. Much attention has been 
focused on the role of entanglement in characterizing 
QPTs in the last few years (e.g. Refs. iSill and ref- 
erences therein). More recently, another approach to 
identify QPTs based on the ground-state fidelity has 
been proposed^ and applied to various many-body sys- 
tems 4^^.1102.13^.1506^7^^9^ Because the fidelity is 

a measure of similarity between states, one anticipates 
that the fidelity should drop abruptly at critical points, 
as a consequence of the dramatic changes in the structure 
of the ground states, regardless of what type of internal 
order is present in quantum many-body states. A per- 
haps more effective indicator is given by the singularity in 
the second derivative of the fidelity (or the so-called "fi- 
delity susceptibility" )^4^^^i^^i^^^^ The main 
advantage of this approach lies in the fact that, since 
the fidelity is a purely Hilbert-space geometrical quan- 
tity, no a priori knowledge of the structure (order pa- 
rameter, correlations driving the QPTs, topology, etc.) 
of the considered system is required for its use. The fi- 
delity approach has been examined in various systems, 
including Dicke model^ one-dimensional XY model in 
a transverse field^^i^ general quadratic fermionic Hamil- 
tonians,®'^ and Bose-Hubbard modelfii^ The success in 
analyzing QPTs in these models shows the generality of 
this procedure. The ground state fidelity is usually dif- 
ficult to calculate, due to the lack of knowledge of the 
exact ground state wavef unctions. Therefore, investiga- 
tions so far are restricted to some particular many-body 
models. Conceivably, in order to understand its validity 
and limitation, a general connection between the fidelity 
and QPTs is highly desired. 

In this paper, we discuss, in a general framework, how 



the fidelity can be related to a QPT characterized by 
nonanalyticities in the derivative of the ground state en- 
ergy. It is found that, under certain conditions men- 
tioned below, the fidelity susceptibility is indeed an ef- 
fective tool in detecting the critical points of first-order 
QPTs (IQPTs) and second-order QPTs (2QPTs), as il- 
lustrated before in several concrete models. However, it 
fails to determine the order of the transition, and may not 
detect higher-order QPTs. We stress that, although the 
fidelity susceptibility can not always detect higher-order 
QPTs, it is possible to signal the Beresinskii-Kosterlitz- 
Thouless (BKT) transition^^ in one-dimensional many- 
body systems, which is a QPT of infinite order. To 
demonstrate this, analytic formulas of the fidelity and the 
fidelity susceptibility are derived for single-component 
Luttinger model, which describes a large class of one- 
dimensional systems possessing a gapless spectrum. The 
fidelity susceptibility is shown to be finite in the thermo- 
dynamic limit even for these gapless systems. Besides, it 
is a smooth function of the controlling parameter driving 
QPTs [see Eq. ([T^ below], as long as the system lying in 
the critical Luttinger-liquid phase. By using these formu- 
las, the fidelity and the fidelity susceptibility can be eas- 
ily calculated for a large class of gapless systems, as long 
as the relation between the Luttinger liquid parameters 
and the controlling parameter driving QPTs is known. 
Employing our general formulas to the one-dimensional 
spin-1/2 XX Z model, we show that the BKT transi- 
tion therein can indeed be signaled by the singularity in 
the fidelity susceptibility. Though we restrict our atten- 
tion to the XX Z spin chain, according to the general 
expressions of the fidelity and the fidelity susceptibility, 
one can expect that same results should apply to all the 
BKT-type QPTs of one-dimensional models. 

We now begin to show the general relation between 
quantum phase transitions and the fidelity susceptibil- 
ity. Only QPTs characterized by nonanalytic behavior 
in the derivatives of the ground state energy are consid- 
ered here. According to the conventional classification, a 
IQPT is characterized by a finite discontinuity in the first 
derivative of the ground state energy. Similarly, a 2QPT 
is characterized by a finite discontinuity, or divergence, 
in the second derivative of the ground state energy, as- 
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suming the first derivative is continuous. Following the 
notations in Ref. 0, the most general Hamiltonian of N 
distinguishable particles governed by up to two-body in- 
teractions can be written as 

^ = E ^^/3 I"') (^^1 + E ^^PyS \^^) 1/3.) (7.1 (^,1 , 



iaf3 



where {jai)} is a basis for the local Hilbert space of par- 
ticle i. For these systems, it has been shown that the 
derivatives of the ground state energy per particle £ can 
be written as^ 



d£_ 
dX 
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tr 
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(1) 



where A is a controlling parameter of the system's Hamil- 
tonian and tr denotes the trace over the degrees of free- 
dom for two particles. U(ij) denotes a matrix with ma- 
trix elements Ua0,^s{ii) = e^^J^^/iV,; -I- V^^^^j, where 

is the Kronecker symbol on particle j, and Ni is 
the number of particles that particle i interacts with. 
pij _ rpj.^_ po{X) is the two-particle reduced density ma- 
trix, which is obtained by taking a partial trace Trj-- over 
all degrees of freedom except particles i and j. Here 
po(A) = l^ol-^)) (^ol-^)! is the density matrix of the 
ground state with j^ol-^)) being the normalized ground 
state. U(ij) includes all the single and two-body terms 
of the Hamiltonian associated with particles i and j. We 
assume that U(ij) is a smooth function of the Hamilto- 
nian parameter A. From Eqs. ^ and one finds that 
the origin of a IQPT (2QPT) is due to the fact that one 
or more of the matrix elements of p^^ {d\p^^) are discon- 
tinuous/divergent at the transition point A = Aci^ 

The quantum fidelity (or the modulus of the overlap) 
F of two normalized ground states |^o(A + 6X)) and 
|4'o(A)) corresponding to neighboring Hamiltonian pa- 
rameters is givenby i^(A-i-JA, A) = | (*o(A-f (5A)|*o(A)) j."^ 
To detect QPTs, a more effective indicator is provided 
by the peak in the "density" of the second derivative 
of the fidelity (or the so-called "fidelity susceptibility") 
5(A)^i^ia^iii^'ii^i^^^ which is free from the arbi- 
trariness of small parameter ^A. In the thermodynamic 
limit (i.e., both the particle number N and the number 
of lattice sites L approach to infinity, while N/L keeps 
finite), 5(A) can be written as^^ 



5(A) = lim lim 



-2\nF{X + dX,X) 
LSX^ 



(3) 



To make comparison with the expressions of the deriva- 
tives of the ground state energy, we rewrite the quantum 
fidelity and the fidelity susceptibility in terms of the den- 
sity matrix po{X) of the ground state. Notice that 



Now expanding po{X + 6X) in powers of SX, and us- 
ing the identity 9>[(^'o(A)|\E'o(A))] — which implies 
Ti[po{X)dxpo{X)] = 0, one can easily show that 



F{X + SX,X) ~ 1 



SX^ 



-Tr 



dpojX) dpojX) 
dX dX 



and thus 



5(A) = lim — Tr 



5/90 (A) dpo{X) 



dX dX 



(5) 



(6) 



F{X + SX, Xf = Tr[po(A)po(A + SX)] 



(4) 



As mentioned before, IQPTs (2QPTs) must come from 
discontinuity in (discontinuity in or divergence of) one 
or more matrix elements of p^^ (dxp^^). Since the ma- 
trix elements of the reduced density matrix p^^ are linear 
functions of those of po, IQPTs and 2QPTs will be associ- 
ated with nonanalyticity in the matrix elements of dxpo, 
and therefore imply the singularity in S{X). That is, the 
singular behavior of the fidelity susceptibility is able to 
signal IQPTs and 2QPTs. This explains the success of 
the fidelity approach discovered in the previous investiga- 
tions. Note that the above conclusion is valid only if the 
discontinuous/divergent quantities do not accidentally all 
vanish or cancel with other terms in Eqs. 0) and (0) [i.e., 
assumptions (b) and (c) in Ref. 1^. However, some lim- 
itations of this approach are discussed in order. First, 
even though 5(A) can be a good indicator of IQPTs and 
2QPTs, it fails to distinguish between them, in contrast 
to the entanglement measurements discussed in Refs. Hi. 
Second, from Eq. we find that 5(A) can not detect 
the higher-order QPTs resulting from the nonanalytici- 
ties in d^po and even higher-order derivatives. Neverthe- 
less, it does not mean that 5(A) always fails to signal the 
higher-order QPTs. Reminding that the two-particle re- 
duced density matrix p'^ is calculated by taking a partial 
trace of po, it is thus possible that the nonanalyticities in 
the matrix elements of d\po cancel one another in obtain- 
ing dxp^K That is, while the discontinuity/divergence in 
p^^ and dxp^-' does imply the nonanalyticities in dxPo, 
the reverse is not true. In this case, dj^E can be contin- 
uous even though 5(A) is singular, and therefore such a 
higher-order QPT can indeed be detected by 5(A). One 
such example is provided by the BKT transition in the 
one-dimensional spin-1/2 XXZ model, where a critical 
anisotropy separates a gapless phase from a gapful phase. 
As demonstrated in Ref. 24, the ground state energy and 
all of its derivatives with respect to the anisotropy are 
continuous at the critical point. That is, it is a QPT 
of infinite order. However, as discussed below, 5(A) does 
become singular despite the regularity of the ground state 
energy at this critical point. 

It is known that many one-dimensional gapless sys- 
tems, which may undergo the BKT transition, can be de- 
scribed by a single-component Luttinger model. They 
include the spin-1/2 XXZ spin chain and the spin-1/2 
J — J' spin chain in their spin fiuid phases, and the 
Bose-Hubbard model in its superfiuid phase, etc. Be- 
fore specifying to the XXZ spin chain, we first derive 



3 



analytic formulas of the fidelity and the fidelity suscep- 
tibility for the single-component Luttinger model. By 
using the bosonization method, the Luttinger model can 
be described by the following effective Hamiltonia n^^i^^ 



dx 



A 



(7) 



Here the bosonic phase field operators <i> and 11 are given 
by22 



k^O ^ ^ 



fc#0 



where Ofc and a^, obey canonical boson commutation re- 
lations: [afc,aj^,] — 5k,k' ■ We note that the Luttinger 
liquid parameters K and u are functions of the control- 
ling parameter A of the original Hamiltonian. When 
K — 1, the effective Hamiltonian reduces to the free- 
boson model, whose ground state \^o{K = 1)) is noth- 
ing but the vacuum |0) of the canonical bosons satisfying 
OfclO) = for all k ^ Q. When X 7^ 1, the effective Hamil- 
tonian contains the unpleasant terms akd-k and a^a^^,. 
These terms can be diagonalized by a Bogoliubov trans- 
formation on the bosons such that 



a-k 



— cosh 6 a_ 



sinh 6 a 



aj. — sinh a_fc + cosh a[. 



with cosh 6* = {l + K)/2VK and sinh 6* = {1-K)/2VK. 
We note that, in the present case, the parameter 9 is in- 
dependent of the momentum k. The transformed Hamil- 
tonian becomes a frcc-boson model in terms of the new 
set of canonical bosons a^- That is, the effective model in 
Eq. ([7]) can be considered as a "quasi-free" boson model, 
where the ground states \^o{K)) for general values of K 
are the vacuum of ckfc satisfying ak\'^(){K)) = for all 
fc 7^ 0. Thus the (not normalized) ground states become 



k^O 



|0) 



(8) 



With these exact expressions of the ground states, now 
one can calculate the ground-state fidelity. From its 
definition, the ground-state fidelity F of two (not nor- 
malized) ground states |^o(^^')) and \'i/o{K)) becomes 
F{K',K) ^ Z{K', K)/y^Z{K', K')Z{K, K), where 
Z{K',K) = \{^o{K')\^o{K))\. By using the expressions 
of the ground states in Eq. it can be shown that 



ZiK',K) = Y[(l- 

fe#0 



sinh 9' sinh 9 
cosh 9' cosh 9 



(9) 



and therefore a general expression of the fidelity is 
reached, 



F{K',K) = l[ 



1 



cosh(6'' - 61) 




(10) 



Here the prime denotes the corresponding variables tak- 
ing their values at K'. Obviously, F = 1 ii K' — K. 
Generically, 1/ cosh(6'' — 9) < 1. Therefore, for systems 
with a large but finite size L, the fidelity in Eq. ([T0|) 
scales as [cosh(0' — 0)]~^, and decays very fast when K' 
separates from K. That is, the fidelity of two different 
ground states becomes zero in the thermodynamic limit, 
despite the two states being in the same phase. There- 
fore, in this case, it is difficult to signal a precursor of 
the QPT simply by seeking a drop in the fidelity. As 
mentioned before, a more effective indicator is provided 
by the fidelity susceptibility given in Eq. ([3|). From the 
analytic expression of the fidelity in Eq. (fTUl) . it can be 
shown that 



lim 



lnF{K + SK,K) 
L 



and one finally obtains the general formula of the fidelity 
susceptibility 



5(A) = 



1 dKV 
4 \K~dx) 



(12) 



Thus we show that S{X) can be finite in the thermody- 
namic limit even for gapless systems. This result agrees 
with the findings in Ref. [13 based on the scaling argu- 
ments. In the present work, we further provide the an- 
alytic expression of 5(A). From this expression, we find 
that 5(A) becomes singular only when dK/dX diverges. It 
is noted that, in terms of the Luttinger liquid parameter 
K, the general expressions of F and 5(A) in Eqs. (fTUl) 
and (jl2|) apply to all one-dimensional systems in their 
Luttinger liquid phases. The precise location of the sin- 
gularity in 5(A) can be determined, once the exact rela- 
tion between the Luttinger liquid parameter K and the 
controlling parameter A driving QPTs is known. 

Now we consider the case of the XXZ spin chain, 
which can be taken as a special case in the Luttinger liq- 
uid description. The Hamiltonian of the one-dimensional 
spin-1/2 XXZ model is written by 



i=i 



A S j S j ^ 2 ) 



(13) 



Here Sj,Sj, and Sj are the spin-1/2 operators at the 
j-th lattice site. The parameter A denotes an anisotropy 
in the spin-spin interaction. The XXZ spin chain is an 
exactly solvable modeli^ It is known that there is a crit- 
ical point of IQPT at A = —1, which corresponds to 
the isotropic ferromagnetic Heisenberg model. At the 
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isotropic antiferromagnetic point A = 1, a BKT transi- 
tion occurs, which is described by a divergent correlation 
length but without true long-range order. 

After applying the Jordan- Wigner transformation and 
bosonization procedure for the spin-1/2 operators, when 
— 1 < A < 1 , the Hamiltonian in Eq. (fT3|) can be mapped 
to the bosonized effective Hamiltonian in Eq. ([7]) f2S, That 
is, the XXZ spin chain can be considered as a "quasi- 
free" boson model. For general values of A (obeying — 1 < 
A < 1), the Luttinger liquid parameters K and u can be 
obtained by comparison with the Bethe-Ansatz solution. 
They are given hy K = {n/2)/[Tr — arccos(A)] and u = 
Wl - A2/(2arccosA).2? For A = 0, we have K ^ I; 
while 7^ 1 for A 7^ 0. Substituting the above exact 
relation between the Luttinger liquid parameter K and 
the anisotropy parameter A to Eqs. and (|12p . the 
expressions of F(A',A) and 5(A) can be obtained. Here 
we discuss 5(A) only, which becomes 

^''^^ " 4[7r - arccos(A)]2 ' ^^^^ 

Therefore, the fidelity susceptibility 5(A) diverges as 
A —> ±1. That is, the singular behavior in 5(A) is able to 
signal either the IQPT at A = —1 or the BKT transition 
at A = 1. According to our previous analysis, the singu- 
larity in 5(A) at A = 1 indicates that one or more of the 
matrix elements of d\po should be divergent at this BKT 
transition. Since the BKT transition is a QPT of infinite 
order, nonanalyticities in the density matrix of ground 
state pq must accidentally all vanish or cancel with other 



terms, such that the ground state energy and all of its 
derivatives with respect to the anisotropy are continuous 
at this critical point. Therefore, the BKT transition in 
the XXZ spin chain does provide an example, where the 
higher-order QPTs can be detected by the singularity in 
the fidelity susceptibility. 

In summary, according to the proposed general relation 
between QPTs and the fidelity susceptibility, the validity 
and the limitation of the fidelity susceptibility in charac- 
terizing QPTs are discussed. Employing our analytic for- 
mulas of the fidelity and the fidelity susceptibility for the 
one-dimensional Luttinger model, it is shown that the fi- 
delity susceptibility can be finite even for critical systems, 
which agrees with the result in Ref. [ij based on the scal- 
ing analysis. Moreover, while the fidelity susceptibility 
may not detect higher-order QPTs in general, we demon- 
strate that the BKT transition, a QPT of infinite order, 
in the spin-1/2 XXZ chain can indeed be signaled by 
the singularity in the fidelity susceptibility. Though we 
restrict our attention to the XXZ spin chain, we believe 
that same results will apply to all the BKT-type QPTs 
of one-dimensional models, such as the transition from 
spin fluid to dimerized phase in the J — J' model, and 
the superfiuid-insulator transition in the Bose-Hubbard 
model at integer filling, etc. 
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